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ABSTRACT 

The muscle spindle, the feedback element In the m/otatlc (stretch) reflex, 
is a major contributor to muscular control. Therefore, an accurate description 
of behavior of the muscle spindle during active contraction of the muscle, as 
well as during passive stretch, is essential to the understanding of muscle 
control. With this goal In mind, animal experiments were performed In order to 
obtain the data necessary to model the muscle spindle. Spectral density 
functions were used to Identify a linear approximation of the two types of nerve 
endings from the spindle A nodal reference adaptive control system was used 
on a hybrid computer to optimize the anatomically defined lumped parameter esti- 
mate of the spindle. The derived nonlinear model accurately predicts the be- 
havior of the muscle spindle both during active discharge and during Its silent 
period. This model will be used In future research to determine the mechanism 
employed to control muscle movement. 


INTRODUCTION 

This paper is concerned with the derivation of a mathematical model which 
describes the behavior of the mammalian muscle spindle of the cat over Its entire 
operating range. The muscle spindle Is a subsystem of a largar system which Is 
responsible for skeletal muscle control. One control loop within this system Is 
known as the myotatlc reflex arc. 

Analyzing the myotatic reflex Itself can be undertaken only when each of the 
* bsystems and element? of this reflex are well defined. It 1$ the goal of this 
p er to experimentally v>d mathematically anal - one of the subsystems of the 
reflex. The subsystem . j» be considered Is me nuscle spindle. 

The muscle spindle is studied during active contraction as wall as during 
Stretch In order to determine the normal functioning of this essential subsystem 
in the myotatlc reflex. An optimization scheme Is utilized on an anatomical 


representation of the muscle spindle In order to escertaln the importance of 
Its various properties. Because It Is Important to relate each mathematical 
term to real components within the biological alemants themselves, the analysis 
of the muscle spindle Is based on Its anatomy rather then treating It as a 
“spaetral fit” of the experimental data. This study utilizes engineering 
concepts to analyze and predict the normal functioning of the mammal 1 an muscle 
spindle. These concepts lend themselves to the analysis of this type of 
biological system because it can be represented as e basic servo mechanism as 
Merton (19521 did in his wall known presentation. 


ANATOMICAL CONSIDERATIONS 




The mammalian muscle spindle Is e complex sensory organ connactad In 
parallel with skeletal muscle fibers, it detects changes In the length of 
the muscle as well as In other perimeters which are related to length, such 
as tension and velocity. The spindle Is composed of several Intrafusal 
muscle fibers contained within one envelope. The multiple fibers contained 
within this envelope fell Into two distinct categories. Of the two fibers, 
the larger, called the nuclear beg fiber, reveals two different regions. 

Tho polar regions of this nuclear beg fiber consist of strletod contractile 
tissue, while the central region Is composed of a noncon tract lie tissue con- 
taining nuclei. The smaller fibers within the envelope are e uniform distri- 
bution of nuclei and striated regions. These fibers ere celled nuclear chain 
fibers due to their chain- like appearance. 

Twi. types of nerve endings Innervate the striated regions of the spindle. 
These are denoted as the ¥. end t, wter fibers. This ge mm a afferent system 
provides the necessary bias to the muscle spindle. Figure I shows the end 
]L motor fibers Innervating the different types of Intrafusal fibers. The 
afferent discharge from the muscle spindle Is of two distinct typos. The 
first, namely, the fe efferent discharge, has Its endings on both the nuclear 
bag and nuclear chain fibers. These afferent discharges are connected mono- 
synaptlcelly with alpha motoneurons. The second type of discharge, the type II, 
has Its ending only on the nuclear chain fiber end Is connected polysyneptlcel ty 
with the alpha motoneuron. 

A more detailed discussion of the anatomy of the muscle spindle Is avail- 
able In the list of references: Matthews [2], Eldred (3), Grenlt fb) and others. 
Of particular Interest are the proceedings of the first Nobel Symposium. 


DATA AQjJISITION 

In order to model a given systam, It must be fully defined by certain 
Information. Knowledge of Input-output relationships Is repaired to formulate 
e mathematical description of the system being studied. Therefore, modeling 
of the muscle spindle requires the recording of Input to the spindle as welt 
as Its output pulse train. Input to the muscle spindle consists of any 
mechanical change In the length of the ext refusal muscle fibers. This change 
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Is transferred to the Intrafusal fibers of the muscle spindle because the 
latter are anatomically placed In parallel with the extrafusal fibers. The 
output of the muscle spindle s/stem Includes the afferent discharge of both 
type la (primary nerve ending' and type II (secondary nerve ending) fibers. 
The following experimental procedures are utilized to obtain the Input-output 
data necessary for modeling the system. 

A simple schematic diagram of the experimental facility explains the 
arrangement of equipment. In Figure 2 the stimulating electrode and monitor- 
ing electrodes are shown on the appropriate nerve bundles. The *tlmulating 
electrode Is orlven by a constant current unit to provide uniform stimulus 
regardless of changes In resistance between the electrode snd nerve bundle. 
The rate and duration of the stimulus pulses are controlled by e stimulator 
which drives the constant current unit. 

Nerve conduction tests are performed In order to differentiate between 
type la and type It nerve fibers from the spindle. A conduction 

velocity of less than 70 m/sec Is generally a type It ending. Further 
i dent I f I cut ton of the nerve fibers Is made by observfrg the dynamic sensi- 
tivity of the fiber output. The type la fiber shows greater sensitivity to 
velocity of muscle movement than the type II ending. 


EXPERIMENTAL RESULTS 

A majority of previous Investigators have driven the muscle by means of 
a mechanical apparatus or external force. Such experiments tend to keep the 
muscle spindle afferent In the range of continuous discharge by controlling 
the gamna efferent bias and mechanical load. It Is assumed that the only 
region of Interest to the physiological system Is that of active discharge. 
This study differs In that It treats the silent period of operation as well 
as the period of active spindle discharge. This Is accomplished by leaving 
the normal nerve pathways Irtact and causing muscle contraction by stimulation 
of the peripheral nerve. The results of these experiments ere discussed In 
detail in succeeding sections. 


Isotonic Loadli 


Isotonic loading of the medial gastrocnemlus-soleus muscle group Is 
accomplished by attachment of a constant tension spring. The measurement 
apparatus provides a low mass mechanism for measurement of muscle movement. 
After attachment o' the spring, the muscle’s length is allowed to stabilize. 
The spindle afferent is monitored In order to ascertain whether the discharge 
has reached steady state. 


Once It has been determined that the muscle spindle and muscl* are at 
steady state, the peripheral nerve to the muscle group Is stimulated with e 
0.5 millisecond pulse at an amplitude sufficient to cause a maximal twitch of 
the muscle fibers. This stimulus puUe and the post-stimulus spindle efferent 
discharge are recorded along with the positional translation of the muscle for 
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a period In excess of one hundred stimuli In order to Insure e large enough 
sample to provide • statistical basis for date analysis. 

The resultant positional translation end spindle discharge ere shown In 
Figure 3. It can be seen that the spindle afferent dltchar^ ^sM^ or Is 
silent, Immediately following the stimulus to the muscle. This silent period 
Is maintained until the direction of the velocity component It r ^ ri9iS - 
Figure 3 also shows that the offset of the silent period Is 
an apparent bursting or high frequency discharge of the spindle efferent. 

The Instantaneous frequency of the spindle afferent discharges provides 
output Information to higher centers In the biological system. Therefore, It 
Is an Important va. latoie which must be measured. Since each post-stimulus 
discharge Is unique, a post-stimulus histogram of Its Instantanoous frequency 
Is necessary (Figure 4), This figure Illustrates the silent period, Its 
offset, and the Instantaneous frequency of the spindle afferent discharge. 


Qbservetlon of static as wall as dynamic results are Important In under- 
standing the normal function of the muscle spindle. Two types of loads, 

Inertial and Isotonic, are utilized to determine the properties of the spindle. 

These different types of loads affect only the muscle 1 • dynamic response end 
consequently, the spindle's dynamic response, while the static characteristics 
of each experiment remain the same. Therefore, the static characteristics for 
both types of load can be discussed together. 

The elastic properties of the muscle ere determined by measuring the 
amount of elongation of the muscle relative to the amount of load. Figure 5 
Illustrates the results of this measurement for four cats. The four curves 
shown ere displaced In the vertical axis but have similar shapes end appear 
as a family of curvas. Tha vartleal displacement Is a natural result of 
normal dlfferaneas between experimental animals. The Important jharactarls- 
tlcs to ba noted here ere the nonlinear nature of the force-displacement 
curves and tha similarity of tha curves to each othar. Tha results Illustrated 
In Figure 5 ere similar to those of previous Invest I gstors. Experiments of 
Wilkie, 1956; Ritchie and Wilkie, 19561 Zejec, 1966; end others thaw an 
exponential relationship batwaon tha length of tha elastic element and the 
force applied. However, the results both of this study end the above works 
disagree with those of Rosenthal (5) who describes the length- tens Ion 
relationship as one of e linear nature with or without the ryotatlc reflex 
present. The data ha p ras an ted are apparently a result of ine dynamic nature 
under which they were obtained. The nonlinear properties of the elastic 
clement shown In Figure $ ere the result of measurements taken with the m^otetfc 
reflex Intact. 

Tha static discharge of the spindle, Illustrated In Figure 6, has character- 
istic s very similar to thoee of tha length- tans I on curvas graphed In Figure 5. 

The points plotted In Figure 6 represent the mean spindle discharge; the normal 
spindle output varies about a mean frequency. Both graphs appear to be logarithmic 
!n nature due to the characteristics of e sharp Initial slope followed by continual 
reduction of slop# as the tension is Increased. 



Inertia l Loading - Dynamic Character! at lea 

A I though the static characrerlstlcs for both types of loads are the 
same, the dynamic characteristics are different. A mass Is used to provide 
a load with Inertia which more closely approximates the type usually en- 
countered In the normal state. The muscle response and corresponding spindle 
afferent are shown In Figure 7. The muscle response to the stimuli shows a 
'leflnlte damped oscillation which results In a secondary silent period In 
the discharge of the spindle afferent. Figure 8 better Illustrates the time 
sequence and duration of this secondary silent period. A family of post- 
stimulus histograms Is Illustrated In Figure $. Whereas the Initial silent 
per od remains approximately constant In duration, the secondary silent 
P ? r l°? ?! n ' fnfsh *s as the load Increases and Is extinguished by loads In excess 
° 3.* k»ioi. Obs-'vation of the family of displacement curves in Figure 10 

shows the seconder) oscillation of the positional data decreasing as the 
' j ncreased * T Ms measure, namely, the sensitivity to secondary 

oscillations, provides a valuable Index of the sensitivity of the spindle 
afferent to displacement, velocity, and acceleration. Another Important 
phenomenon is the relative constancy of the Initial silent period, particularly 
since the muscle dynamics are grossly affected by Increase In load. 


Dynamic Characteristics 

The spindle afferent discharge end corresponding positional translation 
ere illustrated In Figure II. Of Interest here Is the almost complete lack 
of the eyplce* high frequency discharge noted In previous data at the offset 
of the si lari period. The lurst In the afferent discharge readily apparent 
In the type la afferent Is drastically reduced In the type II ending. Its 
afferent discharge frequency decreeses gradually as the muscle returns to 
steady state. Its range of the discharge from the maximum frequency to the 
steady state frequency Is less than twenty percent, whereas the peak discharge 
rate for the type le unding Is four to five times that of the steady state 
discharge. 

Discussion of Experimental Remits 

The preceding sections presented a view of the data generated by the 
various experiments utilized In this study. This section describes the results 
and parameters derived from these experimental Investigations end discusses 
them with reference to the results of previous Investigators. 


The reduction and cessation of pulse discharge during muscle contraction 
or length decrement can be explained by Interpretation of the physical paramaters 
acting on the spindle. The Initial silent period observed in Figure 9 cannot 
be explalnel simply as a response to the mechanical forces acting on the spindle. 
Observation of the large variation In muscle response, illustrated In Figure 10, 
indicates other factors must be present to account for the uniform duration 
of the silent period. Mechanical factors definitely contribute to the silent 
period observed In Figure 9. However, as the load is Increased, the duration 
of the silent period does not decrease. The probable cause of the uniform 
silent period Is Inhibition of fuslmotor neurones by an antidromic volley In 
adjacent neural roots (Eliaway, 1971; Holmgren and Her ton, 1954), This Inhibition 
of fuslmotor neurones Is shown to have a duration of approximately 50 milliseconds 
whl-h coincides with the minimum duration of the silent period found In the 
experimental data. The ihlbltlon of the fuslmotor neurones causes the relaxation 


the silent period observed In 


of the polar regions of the spindle, resulting In ^ r(WJ ^rvec , 

the «*P« r antel data. Jansen and Rudjord [7] also noted the same phenomenon 
and stated. The duration of the silent period of motoneuron was not appreciably 
changed even by large variations In the Inltal tension.'* 


The silent period seen In the datz presented In this paper demonstrates 
the nonlinear nature of the type la and type II endings. An important parameter 
to be determined Is the sensitivity of the silent period to physical forces. 
Comparison of the velocity curve and the discharge occurrence histograms shows 
a relative offset end onset sensitivity of slightly lets than 0.1 centimeters 
per second for offset and approximately 0.9 centimeters per second for onset. 
Utilizing this sensitivity parameter to determine the load required to extinguish 
the secondary silent period, a value of 2.7 kilos Is determined from the velocity 
curve In Figure & Since only t finite number of loeds Is used In these experi- 
ments, Interpolation Is necessary to check on the validity of this value. The 
histograms displayed In Figure 9 show that the secondary silent period Is 
beginning to be extinguished by afferent discharge for a load of 2.270 kilos 
and Is completely extinguished for S.k kilos. This agrees very well with the 
calculated tension of 2./ kilos. As Lennerstrand states, assessment of velocity 
response during length decrement cannot be determined. However, the relative 
values of onset and offset are important In determination of the type of non- 
linearity In the spindle system. The relative vatues of these sensitivity 
parameters are necessary In the prediction of silent period occurrences and 
durations. 


Isotonic Loading - Type II Ending 

The exact role of the type II ending In the reflex system Is not completely 
understood, it is thought, however, that the type II ending provides Information 
to the gamma system which controls the efferent bias to the muscle spindle. The 
date presenter In this section shows that the secondary ending of the muscle 
spindle has very little dynamic sensitivity. Instead, tl ' secondary ending provide] 

a measure of the static position of the muscle. Due to th i characteristic 

the bias control theory appears to have some validity. However, It Is also* 
known that the type 1 1 ending excites flexors end Inhibits extensors and, therefore, 

plays additional roles In the overatl control system. 


Holmgren and Herton found that an Interjected muscle twitch during steedy 
contraction caused a silent period In the electromyogram of the muscle. Before 
reflex action was developed, the motoneuron discharge was already arrested. 

This cassation of motoneuron activity was attributed to antidromic blocking, 

If the twitch was produced by excitation of motor fibers. Similarly, orthodromic 
stimulation produced a depression In motoneuron activity. In our experiments 
It Is apparent that the silent period onset Is attributable to this blocking 
phenomenon, since the onset of the silent period does not depend on the dynamics 
of the muscle system. Verification of this phenomenon Is also given In the 
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experiments utilizing Inertial loads by observation of the secondary silent 
period produced. Clearly the onset of the secondary silent period obeys the 
muscle system dynamics In that case. 

Grant t and van der Meulen investigated the duration of the silent 
period and the cause of the spindle afferent's firing after the cessation 
of activity. The long pause and short pause responses demonstrated were 
attributed to the large differences In length of the spindles contained 
within the muscles studied. The adaptation time of spindle length (approxi- 
mately 10:1 for largest to smallest) to the dynamics of the muscle Is described 
as the predominate cause of the variability In the spindle silent periods. 

This long pause, short pause phenomenon Is not found In the present research. 

In contrast, the silent period described here Is dominated by the dynamics 
of the muscle itself, with the exceptions already mentioned. 


spindle. The model derived Is then tested against tha experimental data 
already presented In order to determine the authenticity of Its behavior. 

Certain constraints, such as approximation of the analog components to 
the anatomical form of the biological elements, must be maintained. in 
addition, tha modal must be able to perform under the Influence of forcing 
functions representing both stretch end contraction. It Is with these 
constraints In mind that tha following models are proposed and developed. 

Logically, the Input or forcing function to the muscle spindle Is the 
displacement of the muscle In which It Is contained. Therefore, the muscle 
Itself must be modeled In order to simulate the Input In the biological system. 
The model of the muscle Is then used as the Input for the model of the muscle 
spindle. 


A relatively constant silent period of motoneuron firing was reported 
by Jansen and Rudjord during single shocks of maximal Intensity to the muscle 
nerve. Three reflex mechanisms were studied In order to determine the cause 
of the silent period duration. These were: antidromic Inhibition, Golgi 
tendon inhibition, and the pause In la excitatory afferent Impulses. It was 
concluded that, under the experimental conditions utilized, the most Important 
factor determining the concoml tant silence of the electromyogram wes the 
pause in la afferent activity.. This agrees with the results shown In the 
presen’ study as the constancy of the la afferent silent period agrees with tha 
reotone* ron silent period. However, this research also defines the overall 
ca^ise <f the la afferent and motoneuron efferent silent periods, since the 
motonei ron efferent Is dependent on the la afferent. 

0% termination of tha frequency of afferent discharge relative to the 
velocity of the muscle is found by measuring the peak velocity of the 
pruparrt Ion and the corresponding peak In afferent discharge. The data for 
thu isc tonic and Inertial load experiments reveal a nonlinear relationship 
between velocity and afferent discharge frequency. The probable cause of 
the differences between the two curves (other than those normally expected 
between experiments) Is the much larger negative velocity of the muscle at 
the onset of active muscle contraction. This Initial large negative velocity 
component apparently overdrives the system. The nonlinear curves Indicate 
that the muscle spindle can respond only up to certain limits, after which the 
gain of the system rapidly decreases. This Is probably due to depolarization 
of the nerve terminals In the spindle. Matthews and Stein also observed this 
nonlinear response in 1969 during sinusoidal changes In muscle length. Their 
experiments were accomplished through the use of mechanical oscillation of 
the muscle, however, and cannot be directly compared to this Investigation. 


THE M00EL 

The previous section briefly presented the results of the experiments and 
analysis of the experimental data. This sactlon utilizes these results In 
determine! Ion of an appropriate model which describes the behavior of the muscle 


Analysts of the spectral density function for the muscle shows that a 
differential equation of at least second order Is necessary to describe the 
response of the muscle. This Is determined from the Increase In spectral 
content prior to the decline of the spectral dans tty function. This behavior 
denotes a complex pair of roots which Implies e second order differential 
equation. Taking this into consideration, the following equation describes 
such a system. 


ax + bx ♦ cx 


f(t> 


. Using mechanical analysis, It Is apparent that the three elements In 
the differential equation represent mass, elasticity, end viscosity. The 
rums Is the lumped sum of the muscle mass end the load to the muscle. 
Elasticity Is the elastic property of the muscle and Its Interconnecting 
fibers, and viscosity refers to the viscous properties of the same elements. 
A model Is formulated from these parameters. 


Elastic Elements 

A reasonable place to begin the modeling process is with the elastic 
elements of the system. Since tha elastic properties of the muscle exhibit 
their characteristics statically as well as dynamically, tha Identification 
schema can be either static or dynamic. Since static characteristics are 
described In tha previous taction and evaluation of tha static properties 
of the element eliminates other dynemlc characteristics of the muscle, static 
evaluation of the muscle fiber Is used. Elasticity Is analogous to a spring; 
therefore. It is a simple matter to model the elastic properties with e 
spring which has the seme characteristic behavior as the experlmr^al date. 
Observation of the data, however, Indicates that a linear spring •« Insufficient 
to provide an accurate description of the muscle’s elastic properties. There- 
fore, several nonlinear equations are evaluated as to their goodness of fit 
to the experimental date. The choice of the function for evaluation Is 
somewhat arbitrary, and perhaps other selections would suffice. However, the 
functions described appear reasonable when observed In relation to the experi- 
mental data. 



linear regression techniques are utilized to determine the optimal 
values of the parameters. Three functions are tested as possible represente* 
tlve lumped parameter estimates of the elastic properties of the muscle 
fibers. The three choices are as follows: 

STRAIGHT LINE ESTIMATE 

y - ax + b straight line estimate 

Q ■ 2^ [VI - (axl + b)J the cost function 

To minimize Q with respect to a and b, the partial derivatives of Q with 
respect to a and b are equated to zero. 

I) ■ 2 Z [Yl - (ax. + b)]X,- 0 

l-l 1 


2) ^ * * 2 2j {Yl - (ax J ♦ b)J - 0 


from these equations are obtained 
n n 

a - 2 X1YI - 7 2 XI 

1-1 1-1 



l-l I-l 


optimal parameter estimate 


ROWER ESTIMATE 


• axb power estimate 


In a ♦ bln x 


y' • a* ♦ bx' by using this transformation the 
previous estimate can be used to determine the transformed parameters. 



In y - In a ♦ bx 

y' - a' + bx the same technique Is utilized to determine 
the transformed parameters. 


The optimal estimates representing the elastic properties of the muscle 
and Its interconnecting tissues are shown in comparison with experimental 
data In Figure 12* Observation of the graphs shows that the exponential 
estimate represents the best fit of the experimental data. The power estimate 
yields a good fit and can also be used as an estimate of the lumped equivalent 
of the elastic elements. The straight line estimate does not represent a 
close enough approximation for use In later work. 

The estimates derived here are used later In the modeling of the muscle 
spindle, since it Is assumed that the striated regions of the spindle possess 
properties similar to those exhibited by the muscle itself. The optimal 
estimates determined are as follows: 

F - 1523.8 X ♦ 0.27 straight line estimate 

F - I0kl.3 X 2 power estimate 

1 .9IX 

F - 168,25 (e - I) exponential estimate 

One term of the differential equation has been determined and results 
In a non: Inear differential equation. The second term to be determined is 
the viscous element In the differentia) equation. 

Viscous Elements 

After determination of the elastic properties of the muscle, the next 
element to be determined represents the viscous property. Viscosity can be 
represented mechanically by a dash pot. Since the viscous properties of the 
muscle can only be exhibited during muscle dynamics, the parameter optimiza- 
tion scheme used must be dynamic In nature. 

The method used to analyze the viscous properties of the muscle is the 
gradient search technique. By Implementing the differential equation of 
the system described earlier on the analog portion of the hybrid computer, 
a dynamic search can be accomplished. Since the properties of the elastic 
element have been defined, one term of the differential equation Is fixed. 

The term to be optimized Is p in the following equation: 

otx + (3 k ♦ cx 2 - f(t) 

To accomplish this end, experimental data Is compared with the data 
generated by the mathematical model on the analog computer. By summing the 
square of the difference between the experimental date and the generated date, 
the error squared cost function Is minimized so that the parameter beta (6) 

Is optimized. The cost function Is defined as: 

J - 2 (XE | - XH,) 2 where XE - experimental data 

l-l 


XM - model data 
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Si nee the optimal -solution Is the point where the cost function (J) 
equals 2 «ro, the best estimate of 8 occurs when the cost Is minimi led. 

The value of the model generated data (XH.) Is the voltage sampled by 
the analog to digital converter at a sample rate determined by the real 
time clock and the nunber of clock count* between samples. The sampling rate 
Is a critical factor. If the sample rate of the digital computer deviates 
too far below the required rate necessary to accurately reconstruct the 
sampled waveform, the minimized value of the cost function would never approach 
zero. Taking these factors Into account, the following algorithm Is implemented 
In the digital computer: 

£ 1 ♦ 1 - $ 1 ♦ J1 

Using the above algorithm, a single parameter search is performed In 
order to determine the optimal value of beta ( 6 ). The resuits of the optimal 
search are illustrated In Figure 1 3 . An observation of the graphs shows that 
an excellent approximation to the experimental data is determined by the 
following nonlinear differential equation: 

mx + 8 * ♦ cx^ * f(t) where m • 1.0 

9 - 0.093 
c • 0.915 

A closer approx I mat Ion can be obtained by replacing the elastic term with 
the exponential estimate previously determined . However, due to considerations 
of ease of implementation on the analog computer, the p over estimate Is used. 

The resultant accuracy of the differential equation In relation to the 
experimental data justifies the deviation from the more optimal estimate of 
the elastic element in the system. 

Husde Spindle Nodal 

One of the original constraints of the modeling approach previously 
stated Is adherence to anatomical relationships within the structure of the 
model . Since some of the elements of the model are already deft pad, the 
relationship of these elements Is determined by referring to the anatomical 
representation of the muscle spindle shown In Figure I. 

Properties similar to those of striated fibers are assumed and, therefore, 
the elastic elements are represented by the same nonlinear spring previously 
derived. The viscous elements are also assumed to possess similar properties 
to those derived for the muscle In which the spindle Is contained. The non- 
contract He portions of the spindle have been shown to have elastic constants 
three tines greatej than that of the striated regions and, therefore, the 
coefficients of those parameters reflect this difference. 

The force generators represented In the lumped parameter model of the 
spindle are assmed to behave in a manner similar to that already demon- 
strated by the muscle containing the spindle. However, since the force 


generators ft(t) generally provide a bias to the system, a steady state 
level Is used to simulate the force of the striated regions of the spindle. 

By modification of the model, this force generation system can Interact with 
the modat in a manner similar to the muscle Itself. 

Secondary Ending 

The tumped parameter model of the nuclear chain fiber of the muscle 
spindle Is shown In Figure 14. The behavior, or Input-output relationships, 
should correspond to the results elicited by the experimental data. In order 
to verify tnese results a comparison of the experimental and mathematical 
data Is a necessity. 

The comparison of experimental results with those obtained from the 
mathematical model Is a better test of the model's validity, since the 
correlation techniques yield only e linear approximation of the nonlinear 
behavior of the secondary ending. The Instantaneous frequency histograms 
and comparison of actual and artificial unit potentials provide sufficient 
basis for verification of the model performance. 

instantaneous frequency histograms of experimental and artificial data 
are shown In Figure 17. The characteristic behavior of the secondiry ending 
is repeated quite accurately by the model for the range of load utilized In 
the experiments. The artificial data Is much more predictable as the random 
variable of fluctuating V efferent Is not present In the model. This 
parameter can be added to the model for future studies. 

Figure 15 Illustrates the similarity of the unit action potential 
post stimulus firing patterns for the actual end experimental situation. 

Tha deviation of the firing pattern of the action potentials of the model 
Is very snail In comparison with that of the experimental data. The reasons 
for this difference are the same as those discussed earlier. The actual 
sensitivity of the secondary ending to muscle dynamics and to K efferent 
bias Is easily adjustable In the model and Is one of the features which 
makes the model useful In studying muscle spindle behavior under a variety 
of system constraints. 

Primary Ending 

Figure presents the tumped parameter model of the nuclear beg fiber 
of the muscle spindle. The following differential equations are derived 
from the model Illustrated In Figure 16. 

•) C|(X - X,) - C,(X, - X 2 ) 2 . B(X| - x 2 > 
b) C|(x, - X 2 ) 2 ♦ B,(X| - X 2 ) - 3C|X 2 2 ♦ B|X 2 



ww solving for X^ fron equation (a) yields 

b,x, - e,x 2 ♦ c,(x - x,) - c,(x, - x 2 ) 2 
X, - x 2 ♦ C, (X - X,) - C, (X| - x 2 ) 2 

57 57 

equation b yields X^ 

2B,X 2 - C | (X | - X 2 ) 2 ♦ B | X , - 3C|X 2 2 

x 2 * J. (X, - X 2 ) 2 a 1 - 1 X, 2 

29 f 2 29) 

The third term shown represents a un Indirect Iona I term In the transfer 
function in as much as the elements represent no load to the system, but 
are directly responsible for the derivative term In the transfer function 
of the type la ending, 

(c) KX3 - B,(X 2 - x 3 ) 


B|X 3 - 9,X 2 - KX 3 



Instantaneous frequency histograms of the model response and experimental 
response are shown In Figure 17. As can be seen, the behavior of the model 
closely follows that of the actual biological system. As In the case of the 
type l! ending, the variability of the post stimulus discharge is markedly 
reduced due to the elimination of the oseudo random variations in the 
discharge frequency of the biological system. 

The post stimulus occurrence histograms shown In Figure 18 also Illustrate 
the ability of the model to approximate the behavior of the muscle spindle 
primary ending. The sensitivity of the model to variations In velocity ts 
adjustable and, therefore, al lows the mods! to depict wide ranges of conditions 
normally found in the physical system. This ability of the model to permit 
adjustment of parameters at lows simulation of a great variety of experimental 
conditions and restrictions. 


The Silent Period 

One of the goals of this study Is to investigate the silent period and 
to predlft Its occurrence and behavior. With this purpose In mind, the 
following silent periods shown In Figure 19 are predicted by the mathematical 
model. Comparison of the predicted silent periods In Figure 19 with actual 
silent periods In Figure 20 Illustrates the accuracy of the projected 
calculations. 

The explanation for the behavior of tha silent period illustrated In 
Figure 19 Is that, once the spindle activity Is extinguished. It will remain 
In this state until It Is forced to resume firing. The cause of the 
resumption in spindle activity shown In Illustration Is a positive velocity 
vector greater then 0.9 centimeters per second. This explains why the dura- 
tion of the silent period Increases rather than decreases with load. 

Initially, ft was thought that the silent period would decrease In duration 
with load until It extinguished due to a reduction In muscle dynamics. However, 
because of antidromic stimulation, the spindle is initially caused to be silent, 
regardless of load, spindle will remain silent for a considerable time 
(approximately 200 milliseconds) unless something causes It to resume firing. 

As load Is Increased, the positive velocity vector exceeds the threshold level 
with Increasing latency. This shift in the time at which the velocity vector 
exceeds threshold Is the ceuse of the longer duration of the silent period 
for increased load. 

T he secondary silent periods for the Inertial load experiments Illustrated 
In Figure 9 do not exhibit the phenomenon of Increasing duration described 
in the previous paragraph. The reason for this deviation In the spindle 
discharge pattern is that the onset of the secondary silent period Is due only 
to the dynamics of the muscle and Is not Initiated by antidromic stimulation. 

The model Is able to predict these silent periods as well, since both onset 
and offset of the spindle activity follow the dynamics of the system. 


IMPLICATIONS OF THE MODEL 

The mathematical model pr *« anted, which was derived fro m experimental data, 
possesses several capabilities not Incorporated In previously published models. 
Consideration of the entire range of muscle functioning el lews the spindle 
afferent discharge to be predicted both during stretch and active contraction. 
Additionally, since the elements comprising the model adhere to the anatomy 
of the muscle spindle. Individual parameter contributions to output of the 
muscle spindle can be studied. Also, sensitivity of the silent period Is 
evaluated and permits the study of the silent period in relation to control 
of muscle movement. 

Since both the muscle and the muscle spindle have been modeled, the 
control law Which governs muscle activity can now be Investigated. We can now 
study the effects of various parameters, such as gamma efferent bias on reflex 
activity. 
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Her ton's follow-up length servo hypothesis [1], which has been widely 
accepted as an explanation of muscle control, does not account for phasic 
muscle movements. His servomechanism can be an adequate mechanism fo- 
postural regulation. 

Considering the information described In previous sections and by 
other Investigators, spindle discharge mediated by % activity Is Inadequate 
to control muscle movement. In addition, the delays of the ¥ route due to 
the propagation time of neural transmission are obvious. A follow-up servo 
would require a high gain In order to maintain Its steady state length, 
regardless of load. However, the ease with which limb posl* ’on is varied 
by application of load verifies the fact that the servo ms chan Ism which 
controls reflex activity has Irw gain. 

A concept which Is taking the place of the follow-up servo theory of 
muscle control is that of servo-essl stance proposed by Hat thews {2]. The servo- 
assistance '.cheme suggests a system which Initiates movement by co-act I vat Ion 
of the cK, jt\4 V pathways. This provides a type of ‘'power steering 11 whereby 
con tree Mon of muscle Is assisted by increased spindle discharge. However, 
Ancep'.ence of this more plausible explanation of the mechanism of motor 
control still leaves many factors obscure. 

How does the CHS process these and the other sensory Inputs such as the 
Golgi tendon o.-gan and joint receptors? The overall control mechanism 
appears to utilize all of these inputs to formulate the control strategem. 

Given adequate estimators of all of the peripheral recaptors and 
actuators, investigation of the role of the central nervous system plays 
In tonic and phasic control of muscle movement can be Investigated. 
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